Influence of surface roughness on the optical properties of plasmonic nanoparticles 
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For plasmonic nanoparticles, we investigate the influence of surface roughness inherent to top- 
down fabrication on the optical properties, and find that it has a surprisingly small influence on the 
position and width of the plasmon peaks. Our experimental observation is supported by simulations 
based on the boundary element method approach. Using a perturbation approach, suitable for 
metallic nanoparticles with a moderate degree of surface roughness, we demonstrate that the reason 
for this lies in motional narrowing where the plasmon averages over the random height fluctuations. 
Surface roughness in large arrays of identical nanoparticles, such as encountered in the context of 
metamaterials, is thus expected to not constitute a major roadblock. 
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Plasmonics bridges the gap between the micrometer 
length scale of light and the length scale of nanostruc- 
turesPJ This is achieved by binding light to coherent 
charge density oscillations of metallic nanostructures, so- 
called surface plasmons, which allow to focus electromag- 
netic radiation down to spots with spatial dimensions 
of the order of a few nanometers P Coupling of quan- 
tum emitters, such as quantum dots or molecules, with 
plasmonic nanostructures can strongly modify their ex- 
citation and emission properties, observable in fluores- 
cence^ or surface enhanced Raman scattering,^ and of- 
fers a unique means for tailoring light-matter interaction 
at the nanoscale. This has found widespread applica- 
tions ranging from (bio)sensoriP and solar cells^ to opti- 
cal and quantum communication technology.^ Plasmonic 
nanoparticles are also at the heart of the emerging fields 
of metamaterials and optical cloakingP 

Huge advances in fabrication techniques over the last 
years allow nowadays to fabricate plasmonic nanostruc- 
tures with well understood and predictable properties. 
Nevertheless, practically all metallic nanoparticles suffer 
fro m size inhomogeneities and nanoscale surface rough- 
ness,* 9 -^ which results in deviations of the plasmonic 
properties from those of idealized nanoparticles A 1 Partic- 
ularly top-down approaches for nanoparticle fabrication 
often involve vacuum deposition of the metal structures, 
which leads to polycrystallinc particles with an apparent 
surface roughnessP^l Despite the important role of sur- 
face roughness, there is still little understanding about 
the impact of such imperfection on the optical properties. 
Recent publications report on the control of nanoscale 
roughness and its strong effect on the nanoparticles far- 
and near-field optical properties! 9 * 10 * 12 ! However, varying 
surface roughness is often accompanied by varying crys- 
tallinity, and therefore the results do not allow for a clear 
distinction between the contributions of the surface and 
the bulk. 

In this paper we provide evidence from experiment, 
theory, and simulation that a moderate amount of sur- 
face roughness has no significant impact on the far-field 
optical properties of metallic nanoparticles. We interpret 
this as a kind of motional narrowing, where the surface 



plasmon averages over the random height fluctuations of 
the metal surface, which leads to destructive interference 
and an overall small net effect. Our findings might be 
beneficial for the design of metamaterials, which rely on 
large ensembles of practically identical particles, as well 
as devices based on the far-field properties of plasmonic 
nanoparticles. 

Experiment. — We measured the scattering spectra of 
individual poly crystalline (rough) gold nanoparticles to 
probe the influence of unavoidable surface roughness on 
the plasmonic signature of nanoparticles of nominally 
identical shape. The particles were fabricated by electron 
beam lithography on an indium-tin oxyde (ITO) covered 
glass substrate. The substrate was coated with a polymer 
resist, which was then exposed and chemically developed. 
Vacuum deposition of gold and a liftoff process lead to 
polycrystalline particles (crystallite size ~ 20 nm) of de- 
signed shapes,^ but with apparent surface roughness, as 
shown in Fig. |T^c) . The particles were probed in a dark 
field microscope, collecting the scattered light by a 40 x, 
0.75 numerical aperture objective and analyzing it with 
a spectrograph. 

Fig. [ija) shows exemplarily scattering spectra of four 
individual 50 nm wide, 100 nm long, and 45 nm high gold 
nanorods. Although the exact particle shapes differ due 
to nanoscale roughness [Fig. [l[c)], the dipolar plasmon 
resonance positions and the spectral widths of the four 
particles under consideration differ less than 20% of the 
full width at half maximum. Similar results were also ob- 
served for other nanoparticles shapes (not shown). From 
this we conclude that surface roughness has a surpris- 
ingly small influence on the optical properties of plas- 
monic nanoparticles, at least in the regime where the 
roughness does not lead to a noticeable change of the 
particle aspect ratio. 

Simulation. — We additionally performed simulations 
based on the boundary element method (BEM) j^HUl us _ 
ing a dielectric function representative for golop31 and an 
effective refractive index of rij, — 1.65 for the ITO covered 
glass substrate. The particle shapes in the (a;, y) -plane 
were extracted from the SEM images [Fig. 1(c)]. We 
use a particle height of 45 nm and round off the edges 
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FIG. 1: (color online) (a) Scattering spectra of four individ- 
ual gold nanorods excited with a polarization parallel to the 
particles long axis. Since the repeatability of the scattering 
intensity is only within ~ 20%, the spectra are normalized 
to facilitate the spectral comparison, (b) Spectra simulated 
within our BEM approach, using the gold dielectric function 
of Ref. [H] and a refractive index of rib = 1-65 for the sub- 
strate, for the particle shapes extraced from the SEM images 
reported in panel (c). The particle height is 45 nm. The 
lengths of the scalebars in the lower panels are 100 nm. 



with a curvature radius of 5 nm. The simulated scatter- 
ing spectra are reported in Fig. 1(b). In accordance to 
experiment, we find a surprisingly small influence of the 
roughness-related particle shape on the plasmon peak po- 
sitions. Also the absolute peak heights vary by less than 
5 percent. We note that there are small differences be- 
tween experiment and theory, such as the weak shoulder 
around 2 eV, but the overall agreement regarding both 
peak positions and widths is striking. 

The agreement between the measured and calculated 
widths of the plasmon peaks is quite remarkable since 
in our simulations we do not consider plasmon dephas- 
ing due to surface roughness scatterings, but only in- 
clude electrodynamic decay channels. From supplemen- 
tary calculations we find that the absorption cross sec- 
tions are about three times larger than the scattering 
cross sections. Thus, ohmic losses (described through 
the imaginary part of the dielectric function) dominate 
over radiative damping as well as over dephasing losses 
due to surface roughness scatterings. For smaller parti- 
cle sizes, comparable to the electron mean free path of 
several nanometers, such dephasing will become increas- 



ingly important and could no longer be neglected in the 
simulations. 

To inquire into the reasons for the almost negligible 
influence of surface roughness on the plasmon peak po- 
sitions, in the following we use cylindrical nanorods as a 
showcase system, because of their extreme sensitivity to 
variations of the dielectric environment an d shap ed We 
also employ the quasistatic approximationp^D w hich is 
justified for small particles and will allow us further be- 
low to introduce a perturbation analysis. With respect 
to the retarded BEM simulations, the quasistatic approx- 
imation is expected to lead to a slight overestimation of 
the plasmon peak energy. 

In our simulations we model surface roughness by 
adding stochastic height variations to the smooth sur- 
face of an ideal nanoparticle. In two dimensions and for 
a box with periodic boundary conditions, height varia- 
tions with a Gaussian autocorrelation can be obtained 
by attaching to all Fourier coefficients arbitrary phase 
factors e l( ^ rnd viz. 



h(x,y) = A/i9?e 
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where T~ x denotes the inverse Fourier transform and 
<j\ is the variance of the height fluctuations. We next 
map h(x, y) to the nanoparticle surface and displace the 
vertices of the nanoparticle along the surface normal di- 
rections!^ A typical realization of surface roughness is 
depicted in Fig. [2jb) for a nanorod. 

Figure |2ja) shows spectra for gold nanorods for differ- 
ent height-to-diameter ratios. Throughout we set the rod 
diameter to 30 nm and use ah = 3 nm. For each spectrum 
we average over hundred random realizations of surface 
roughness and assume light excitation with a polarization 
along long axis of the rod. One clearly observes that when 
increasing the amount of surface roughness, by choosing 
different values of Ah, the shape of the averaged spec- 
tra does not change (with exception of an overall small 
red shift). Only the variance of the individual spectra 
increases with increasing Ah, as indicated by the error- 
bars. Thus, our numerical simulations demonstrate that 
surface roughness has a surprisingly small influence on 
the optical properties of plasmonic nanoparticles. 

Things change considerably if we use the more regular 
height variations of a rippled rod, depicted in Fig. [2] 
that were recently suggested as a viable model for surface 
roughness.^ Here the spectra strongly shift to the red 
when Ah is increased, as depicted by the gray lines on the 
left of Fig. [2] Thus, it appears that the irregular shape 
of the surface height fluctuations is responsible for the 
overall small shift of the plasmon peak position. Indeed, 
when ripples are introduced not only along the symmetry 
axis of the nanorod but also around its circumference 
(here 6 ripples) the redshift becomes strongly reduced, 
as shown by the dashed gray lines. 

Theory. — We next develop a perturbation analyis 
which will help us to understand the effects of surface 
roughness more deeply. Let us first recall the basic ele- 
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FIG. 2: (color online) (a) Simulated spectra for nanorods with 
a diameter of 30 nm and for height-to-diameter ratios of 2 : 1, 
3:1, and 4:1. In our simulations we use a dielectric func- 
tion representative for golcP^ and a homogeneous dielectric 
background (nj = 1.5), and add stochastically surface rough- 
ness according to the prescription of Eq. ([I]) with 0-^ = 3 nm 
and for different Ah values, which determine the height vari- 
ations of surface roughness. Each spectrum is averaged over 
100 randomly generated nanoparticles. The solid lines report 
the spectra of the smooth nanorods in the panel for Ah — 2 
nm, and the mean values of the spectra otherwise. The gray 
lines report results for rippled rods, as discussed in the text. 
The spectra for ratios 2 : 1 and 3 : 1 have been magnified 
by factors of 15 and 2, respectively, (b) Surface charge dis- 
tributions a a for the optically active plasmon modes for the 
smooth (left), rippled (middle), and rough (right) nanorods. 

merits of the BEM approach. For a metallic nanoparticle 
with dielectric function e m (oj) embedded in a dielectric 
background with constant e&, the solutions of the Pois- 
son equation within the two media are given by the Green 
function G(r, r') = l/\r — r'\. The electrostatic potential 
can be written in the ad-hoc form (j> — (G, a) + ^> e xtp^ 
where a is a surface charge distribution located at the 
boundary dV of the metallic nanoparticle, which has to 
be chosen such that the boundary conditions of Maxwell's 
equations are fulfilled, (G, a) = f gv G(r, s)a(s) ds de- 
fines an inner product, and 4> ex t is the potential of the 
external perturbation. 

Continuity of the normal component of the dielectric 
displacement at the metal-dielectric interface gives an ex- 
pression 

A w cr + I — ,cr \ = — — , A = 2tt 2 

\ dn I dn e m - e b 

that can be used for the calculation of the surface charge 
distribution. Here dG/dn = F and d<fi cxt /dh = (f>' cxt 
denote the surface derivatives of the Green function and 



the external potential, respectively. 

We next define the right and left eigenvectors a k 
and a k of th e surface derivative of the Green function 
througfPHH 

(F,a k )=X k a k , (a k , F) = X k a k , (3) 

which form a biorthogonal set with (a k , ov ) = dkk' ■ The 
functions a k can be interpreted as the surface plasmon 
eigenmodes, and the response to any external perturba- 
tion can be decomposed into these modes viz. 

— (4) 

Apparently, a given mode k gives a noticeable contribu- 
tion only if the coupling (a k , (j)' cy . t ) to the external poten- 
tial (in our case plane wave illumination) is sufficiently 
strong and if the denominator becomes small. The plas- 
mon resonance condition translates to Re [A(cj) + \ k ] — 
0, when assuming that the imaginary part of e has an 
only weak frequency dependence. This is an extremely 
useful expression as it allows to separate the structural 
properties (described by X k ) from the material properties 
(described by A). 

We are now in the position to analyze the effects of 
a moderate surface roughness, which we model as a dis- 
tortion of the surface dV from its ideal shape. In turn, 
the surface derivative of the Green function F changes 
to F + SF, where SF is expected to have the same ran- 
dom character as the surface fluctuations.^ How does 
X k , which determines the peak positions of the plasmons, 
change in case of surface roughness? When SF is suffi- 
ciently small we can employ perturbation theory, in com- 
plete analogy to quantum mechanics. We treat F as the 
unperturbed part and SF as the "perturbation" . Follow- 
ing Ref. [55] we introduce for a given plasmon mode a the 
projector 

n ,0x0 
V0 _ v-^ cr k a k 

a -Z. A o_ A 0> 1 & J 

k^a a k 

where the superscript indicates the eigenvalues and 
cigenfunctions for the ideal nanoparticle surface. The 
corrections to A a then become within lowest order per- 
turbation theory 

A* = (<7°, SFaty (6a) 

Xl = (a° a ,6F^6Fa° a ). (6b) 

From this we can draw an important conclusion for 
the dipolar modes, whose <r" distributions usually vary 
smoothly on the length scale of surface roughness. Since 
Eq. (|6^,) averages the random variations SF over the un- 
perturbed distributions <r°, the positive and negative 8F 
values associated with the height variations of surface 
roughness will become effectively averaged out. This re- 
sults in a small correction, and consequently surface 
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FIG. 3: (color online) Comparison of full simulations and per- 
turbation theory of Eq. for the positions of the resonance 
energies of the optically active surface plasmon modes with 
polarization along the long axis of the rod. The open and 
filled symbols report results obtained within first and second 
order perturbation theory, respectively. In the inset we show 
a blow-up for the nanorod with a height-to-diameter ratio of 
4 : 1. 



roughness does not affect significantly the optically active 
plasmon modes. 

In a sense, this finding is similar to motional narrowing 
in semiconductor quantum wellsj^ where the motion of 
an optically excited electron-hole pair (exciton) is subject 
to the potential induced by the local monolayer fluctua- 
tion inherent to quantum wells. When the exciton propa- 
gates through the well it "averages" over the fluctuations 
of the random potential, which results in a narrowing of 
the exciton lineshape. 

In Fig. [3] we compare the results of our full simula- 
tions with the predictions of Eq. ^ . The open and full 
symbols show the true resonance energies (obtained from 



Re [A(cj) + Afe] = 0) and the perturbation results (ob- 
tained with A°+A^ and A°+A^ + Aq, respectively). Sym- 
bols on the diagonal correspond to the situation where 
perturbation theory and the full BEM simulations coin- 
cide. In all cases the shift of the resonance with respect to 
the positions of the smooth nanorod (dashed horizontal 
lines) is small. Quite generally, first order perturbation 
theory, which ignores any variations of the surface charge 
distribution (a a = <j°), leads to a blue shift. Only for 
second order perturbation theory, which includes modi- 
fication of a a , the red shift of the plasmonic resonances 
is properly reproduced. 

As for the rippled rod, we observe that the first-order 
corrections A* are small. This is because the argument of 
motional narrowing in principle also applies here. How- 
ever, through the admixture of excited surface plasmon 
modes c , described by the second-order correction A^, 
the surface plasmon can accomodate to the regular height 
variations of the rippled rod, and the plasmon peak po- 
sition becomes strongly red-shifted. No corresponding 
conclusions prevail for the stochastic height variations. 

Summary. — In summary, we have investigated the in- 
fluence of surface roughness on the optical properties of 
plasmonic nanoparticles, and have found a surprisingly 
small effect. Using a simulation and perturbation theory 
approach, we have been able to trace back our findings 
to a motional narrowing, where the plasmon averages 
over the random height fluctuations. As no correspond- 
ing conclusions prevail for the near-field optical proper- 
ties, our results are in accordance with the findings of 
"hot spots" in fluorescence or surface enhanced Raman 
scattering experiments. 

This work has been supported in part by the Austrian 
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